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ABSTRACT 



Aims. This study seeks to provide empirical constraints on the different physical components that can act to yield temporal variability in pre- 
dominantly or partially wind-formed optical lines of luminous OB stars, and thus potentially affect the reliable determination of fundamental 
parameters, including mass-loss rates via clumped winds. 
Methods. Using time-series spectroscopy from epochs spread over ~ 4 years, we present a case study of the 09.5 supergiant a Cam. We 
demonstrate that the Hel /15876 (2'P''-3^D) line is an important diagnostic for photospheric and wind variability in this star The actions of 
large radial velocity shifts (up to ~ 30 km s"') in the photospheric absorption lines can also affect the morphology of the H„line profile, which 
is commonly used for measuring mass-loss rates in massive stars. 

Results. We identify a 0.36-day period in subtle absorption profile changes in Hel A5S76 , which likely betrays photospheric structure, perhaps 
Q ■ due to low-order non-radial pulsations. This signal persist over ~ 2 months, but it is not present 2 years later (November 2004); it is also not 
^ seen in the stellar wind components of the line profiles. Using a pure Ho-line-synthesis code we interpret maximum changes in the red-ward and 
C/2 peak emission of a Cam in terms of mass-loss rate differences in the range ~ 5.1 x 10"* to 6.5 x 10"* M© yr"'. However, the models generally 



Oh 



fail to reproduce the morphology of blueward (possibly absorptive) regions of the profiles. 

Conclusions. The optical line profiles of a Cam are affected by (i) deep-seated fluctuations close to, or at, the photosphere, (ii) atmospheric ve- 
locity gradients, and (iii) large-scale stellar wind structure. This study provides new empirical perspectives on accurate line-synthesis modelling 
of stellar wind signatures in massive luminous stars. 
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1. Introduction 

Through their powerful stellar winds, high ionizing fluxes, and 
ultimate demise as supernovae, massive stars have a profound 
effect on the dynamical and chemical evolution of our own and 
other galaxies. The hot star winds are driven by line scatter- 
ing of radiation, and the extent of the mass-loss is a major 
factor in the evolution of these stars; it can for instance influ- 
ence whether the e nd-state is a neutron star or a black hole (e.g. 
iHeger et a! l2003l) . Therefore it is important that the nature of 
these outflows is well understood, and in particular the mass- 
loss rates are reliably determined. 

One of the most serious obstacles in this path arises 
from detailed optical and ultraviolet observations that re- 



veal widespread line profi 
structure in the winds (see 


e variabihty, implying clumped 


"Prinia. Massa & FuUerton', "20021 


IMarkova et all 12004; 


Morel et al. 




2004| plus references 



within). The developing picture of massive star winds as a 
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clumped medium is furt her supported by rec ent model atmo- 
sphere calculations; e.g. 'Bouret e t alJ ll2005t ) need to include 
substantial clumping factors in their models to match key wind 
features in the spectra of O-type stars, implying significant 
downward revisions in the canonical mass-loss rates. Clumping 
may be a consequence of the inherent instability of the radia- 
tion driving mechanism of the wind, or of the creation of coher- 
ent structures arising from photospheric disturbances, or both. 
In either case, wind clumping is inextricably linked to the fun- 
damental physical problem of radiation driven winds. 

Our objective in the study presented here is to exploit de- 
tailed time-series optical observations of the late-O supergiant 
a Cam (HD 30614) to provide empirical constraints for ab ini- 
tio model atmosphere codes that are commonly used to syn- 
thesise spectra i n order to de rive fundamental stellar and wind 
param eters (e.g lHilliereTSillT 998, CMFGEN.Pauld rach et all 
I2OOIL WMbasic'Puls et al.''2005, FASTWIND). Our approach 
is to constrain the components that contribute to changes evi- 
dent in diagnostic wind and photospheric lines. We seek to 'de- 
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Table 1. Fundamental parameters for a Cam. 



Table 2. Summary log of observation. 



Parameter 



Value 



Reference 



Observatory UT Date 



HJD 2 450 000 + 



No. of 
spectra 



Spectral Type 09.5 la 



Walborn (1973) 



RJRe 
Vg sin(i) 

P,„t(max) 



29000 
32.5 

115kms-' 
20 

1550 km S-' 
-14.3 days 



Repolust et al. (2004) 
Repolust et al. (2004) 
Pennv (1996) 
Schaller et al. (1992) 
Haser (19957 



compose' the variability patterns and assess in particular the 
degree to which fluctuations in the wind-formed features are 
due to 'contamination' by photospheric and deep-seated atmo- 
spheric activity. This is particularly critical in late O and early 
B stars, where the Hq^ emission may be weak; see e.g. the im- 

S irints of photospher ic signals seen in the B aimer lines of e Ori 
Prinia et al. 1120041. BO la,). Our study of a Cam is also de- 
signed to probe causal links between photospheric disturbances 
and spatial wind structure. 

The target a Cam is previously known to exhibit line- 
profile vari a biUty i n H„ (e.g. JEbbets ( 1982); FuUerton ( 1990); 
iKaper et a] ll 19971) ) as well as in some other absorption 
lines in the optical spectrum frullert on. Gies & BoltonLri996 f. 
[Markova, 2(302). Whist its UV resonance lines are strong and 
saturated. iLamers et alJ ( Il988l) reported on weak low- velocity 
fluctuations and changes in the soft blue wings of saturated 
('black') absorption troughs. Some fundamental parameters of 
a Cam are listed in TableQ] 



2. Observations 

Our investigations are based on three datasets; (i) time-series 
of primarily He 1 /15875.67 obtained during January, February 
and March 2002, (ii) time-series of Hq,, Hel /16678 and 
Hel /15 876 obtained during November 2004, and (iii) compar- 
isons of (i) and (ii) to H(,observations secured in 1999 and de- 
scribed by Markova ( 2002). A summary log of the observations 
of a Cam is given in Table |2 

The majority of the 2002 data were obtained using the 
Coude spectrograph of the 2-m telescope at the National 
Astronomical Observatory, NAO (Bulgaria) equipped with a 
PHOTOMETRICS CCD (1024x1024, 24^^). ' The spectro- 
graph was used in two different configurations producing spec- 
tra with a reciprocal dispersion of ~0.2 A pixel"' (effective res- 
olution R=15 000) and of ~0.1 A pixel"' (effective resolution 
R=30 000) over a wavelength range of respectively, -200 and 
~ 100 A . The typical signal-to-noise of the individual spectra, 
adjacent to Hel /15876 , ranges between 100 to 500. 

The 2002 time-series from NAO is complemented by ob- 
servations obtained between 21 to 23 January, 2002 from the 
0.9-m telescope at Catania Observatory (CO, Sicily). The setup 
employed an Echelle spectrograph and a SlTe CCD detector, to 



' This detector was characterized by an rins read-out noise 3.3 elec- 
trons per pixel (2.7 ADU with 1.21 electrons per ADU) 



Catania 0.9m 2002 Jan. 21,22,23 2296.400-2298.577 50 

NAO 2m 2002 Jan. 25,26,27 2300.240-2301.396 64 

NAO 2m 2002 Feb. 25,26 2331.246-2332.486 36 

NAO 2m 2002 March 2,3 2336.220-2337.433 30 

NAO 2m 2004 Nov. 26,27,28 3336.367-3338.654 29 



yield a spectral resolution R ~ 20 000 and a typical signal-to- 
noise of ~ 130. 

The 2004 time-series (Table 2) were also secured at 
the NAO 2-m telescope employing the Coude spectrograph 
to cover wavelength regions including Hel /15876, Hq. and 
Hel /16678 . The effective spectral resolution of these data is 
15 000 with signal-to-noise (adjacent to Hel /15876) between 
300 to 400. 

The data from both observatories were reduced in a ho- 
mogeneous manner, with standard procedures for bias sub- 
traction, flat-fielding and wavelength calibration. The line pro- 
files of a Cam were subsequently normalised by fitting a 
low-order polynomial through selected continuum windows. 
Finally, the telluric water vapour lines were removed by di- 
viding indi vidual spe ctra with a scaled model telluric spectrum 
(e.g. Markova & Valchev. .2000) . Note that the internal stabil- 
ity of the velocity scale in the extracted spectra (specified by 
the standard deviation in velocity of the IS Nal D line at A 
5889.95) is less then 1 km s 'for both the NAO and CO data 
sets. Throughout the paper we have correc ted to the stellar rest 
frame for a radial velocity of 6 km s"' (e.g. lHoffleit & .TaschekL 
fl982) . 

The intensive Hel /15876 time-series of 114 spectra ob- 
tained between 21 to 27 January 2002 is therefore our primary 
dataset for probing short-term modulated or periodic line pro- 
file variability. 

3. Variability characteristics 

With the possible exception of Hq., the optical spectrum of 
a Cam is typical for its late O spectral class, with (mostly) 
symmetrical He I absorption lines, very weak metal lines (e.g. 
C IV /1/15801, 5811), and weak He II ^4686 (absorption). We 
are confident that ours is a case study of a 'normal' late O-type 
supergiant and focus in this section on the temporal behaviour 
of Hel ^5876 and Hq. 

3.1. Hel A5876(2^P^^-3^D) 

This is an important diagnostic line, since the lower level of the 
transition can become substantially populated in the spectra of 
massive stars, which then effectively becomes the ground state 
and the line may reveal P Cygni-like signatures. 

The night-to-night variations are contrasted against 
shorter time-scale (hourly) changes in Fig. 1. Even though 
Hel /15876is predominately a photospheric line in a Cam, the 
night-to-night mean profiles in Fig. 1 (upper-left panel) reveal 
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Fig. 1. Examples of variability in the Hel /15876 and HaUne profiles of a Cam are shown for two characteristic time-scales; (i) 
night-to-night mean profiles (upper panels) and (ii) ~ hourly changes within a single night (lower panels). 
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Fig. 2. The temporal variance spectrum (TVS) for the 
cases of night-to-night changes (solid lines) and ~ hourly 
changes (dashed lines). The horizontal lines in the upper 
panels indicate the respective 95% confidence levels. The 
Hel /15876 fluctuations reward of +500 km s"' are an artifact 
from the removal of telluric lines. The vertical lines (bottom 
panel) mark +Vg sm{i). 



evidence of a variable stellar wind component, with absorptive 
changes in the blue- wing which contrast with the less extended, 
mostly 'stable' red-wing and a weak emission. The data indi- 



cate that the deep-seated inner wind region (< 2 R*) of the star 
is not steady This night-to-night variability is accompanied by 
more rapid and subtle changes (down to ~ 2-3% of the con- 
tinuum level) in the core absorption. This behaviour is seen in 
intensive (hourly) time-series captured during individual nights 
(Fig. 1; lower-left panel) 

To quantify the significance level of the line profile vari- 
ability as a function of velocity, we have determined the 
Te mporal Variance Spectrum (TVS ) according to the methods 
of iFuUerton. Gies & Boltoiil (ll996V The results are shown in 
Fig. 2 (left-hand panel). The night-to-night stellar wind vari- 
ation in Hel /15876 extends over ~ -500 to +200 km s '(i.e. 
~ -0.3 Vco to +0.13 Vco). The ~ hourly core absorption changes 
are essentially symmetric about rest velocity and constrained 
within the boundary marked by + sin(/)thus suggesting the 
presence of photospheric activity in the star. 



3.2. Ha 

The corresponding results for H^are shown in the right-hand 
panels of Figs. 1 and 2. The mean H„profile (e.g. Fig. 2) con- 
sists of a well-developed P Cygni-like core superimposed on 
broad emission wings. In O-type supergiants of same spectral 
type, H„ Profiles with similar morphology to a Cam are in- 
dicative of relatively hig h wind density connected to a higher 
mass-loss rate (e.g. Markova et allbOOSI) . The well developed 
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Hq, profile observed here may (in part) reflect the la luminosity 
class, and perhaps also the presence of clumps in the outflow. 

From Figs 1 and 2 it is also obvious that the emission peak 
of H„ in a Cam is red-shifted, typically occurring at +100 km 
s Red-shifted emission has also been observed in UV reso- 
nance lines of O-type stars, as well as in H„. These two phe- 
nomena, though similar on a first glance, seem to have different 
interpretations. While in the former case the red-shifted emis- 
sion can be explained in terms of "micro-turbulence effects" 
with v„,icro of the order of O.OlvoolHamann, 1980), in the later 
case it likely results from the interaction between the red-wing 
of the Stark-broadened photospheric profile and wind emission 
of cer tain strength, namely log < p > > 13.6^ (Markova et a| 
I2OO5I) . Indeed, with its "mean wind density", log < p >, 
= 13.64, a Cam is just above this limit and thus consistent with 
the interpretation suggested by Markova et al. 

The TVS for Ha nightly means (Fig. 2) extends almost 
symmetrically in velocity space (blue and red limits at about 
+600 km s"', i.e. 0.39voo) with blueward amplitudes (with re- 
spect to the rest wavelength) slightly stronger then the red-ward 
ones. As can be seen from a comparison to Figure 1, the ob- 
served blue-to-red asymmetry in the strength of the TVS am- 
plitudes is more likely due to absorptive blue-ward H^features, 
perhaps caused by wind structure localised along the line-of- 
sight to the stellar disk. In contrast, the overall emission profile 
is produced by gas occupying a greater volume than the gas 
responsible for (episodic) absorption features. 

Generally the H„ profile of a Cam is not substantially vari- 
able on hourly time-scales. However, an exceptional case is 
shown in Fig. 1 (lower-right panel), where localised systematic 
changes are clearly evident over ~ 10 hours. These intriguing 
variations are tightly constrained within +Vg sm(i) (s ee Fig. 2). 
Such a behaviour has previously been documented bv lMarkoval 
ll2002h and we discuss it further in Sect. 3.3.2. 

In a given epoch, over nightly time-scales the total H„ 
equivalent width may vary by up to 10%; however, differences 
of 30-40 % are also noted between line profiles secured in dif- 
ferent years (see Sect. 4). 

3.3. Photospheric activity 
3.3.1 . The 2002 time series 

It is clear from Figs. 1 and 2 that the temporal stellar wind 
activity in a Cam dominates the overall variance levels of 
Hel /15876 and H^. Since the wind changes occur over time- 
scales of days, we are still able to pursue the effects of shorter 
time-scale photospheric variability in these lines (particularly 
in Hel /15876). To emphasis the photospheric changes, the in- 
dividual line profiles were normalised by the corresponding 
nightly mean profile, to 'remove' the larger amplitude night- 
to-night changes, and time-series analyses was then carried out 
on the quotient spectra. 

The most suitable dataset for a search for short-period sys- 
tematic variability is the intensive Hel /15876 time-series ob- 

^ < p > is the "mean wind density"calculated at a typical location 
of 1.4 R, 
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Fig. 3. Two-dimensional representation of Fourier analysis of 
the Hel /I58761ine profile variability. The significant period in 
the power spectrum (right-hand panel) is at ~ 0.36-days. (The 
vertical lines mark +Vj, sin(i).) 

tained from Catania and the Bulgarian NAO between 21 to 
27 January, 2002, totalling 114 spectra (Table 2; note that 
HffWas not observed during this setup). We carried out a peri- 
odogram analysis of these data, using the CLEAN algorithm 
(Roberts et al., 1987) to deconvolve the features of the win- 
dow function from the discrete Fourier transform. The two- 
dimensional grey-scale representation of the periodogram for 
He I is shown in Fig. 3 (for a gain=0.5 and 100 iterations; lower 
values of gain or greater iterations yield essentially the same 
power spectrum). Though the power across the absorption core 
is slightly dis-jointed, the results indicate a (dominant) primary 
frequency of ~ 2.78 days corresponding to a period of 0.36 
days. (N ote this period is consistent with the modulations re- 
ported bv lMarkoval ( I2OO2I) in data for Hel /16678 secured in late 
1998 and early 1999.) We are very confident that 0.36-d period 
is not directly connected to the sampling window of each night, 
where the length of the nightly runs in 2002 is between ~ 0.12 
to 0.27 days (Table 2). Note, that a Fourier analysis of the origi- 
nal He I data (i.e. not residual with respect to the nightly mean) 
reveals essentially the same modulation signal as in Fig. 3, but 
does not show any significant power on this (or any other) pe- 
riod in the more blue-ward region dominated by night-to-night 

wind variability, i.e. over 100 to -500 km s"' . 

Grey-scale images of the phase versus velocity behaviour 
of the individual residual He I spectra on the 0.36-d period are 
shown in Fig. 4. The left-hand panel show the coherent be- 
haviour during our primary (January 2002) time-series. The 
sparser datasets secured in February 2002 and March 2002 
(Table 2) are not suitable for a reliable Fourier analysis, but 
they do provide some indication that systematic modulations 
persist for several weeks over comparable timescales, though 
not with the same phase relationship. 
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Fig. 4. Grey-scale representations of variability in Hel /l5876for individual spectra phased on the 0.36-d period. Individual 

quotient spectra with respect to the mean are shown, such that darker shades indicate regions of enhanced absorption (reduced 
flux) compared to the mean. The total range of intensities plotted is 0.99 (black) to 1.01 (white). 



The data from January 2002 reveal some indication for pro- 
gressive changes across the absorption profile, though the ac- 
celeration of the pseudo-absorption (emission) feature is not 
substantial. There is evidence for characteristic 'blue-to-red' 
motion (see Fig. 5), accommodated within the projected ro- 
tation velocity (115 km s"'). We estimate a prograde feature 
travelling across the line centre with (dV/d(^) ~ 80 km s~ Vcycle 
(Fig. 5). This behaviour would be consistent with a low-order 
sectorial non-radial pulsation mode. These time-series of a 
Cam are however not extensive enough, nor sufiiciently high 
signal-to-noise, to attempt more detailed modelUng to deter- 
mine pulsational parameters. 



The wavelength range of the Bulgarian NAO data of 
January 2002 includes the very weak photospheric metal fines 
of C IV /1/15801, 5812. Unfortunately, the intensity of these 
lines is rather small but we are confident that they are also tem- 
porally active. The C IV /15801 data for 26 January, 2002 cover- 
ing ~ 5.3 hours is shown as a grey-scale ('dynamic') represen- 
tation in Fig. 5, together with the corresponding Hel /15876 line 
profiles. The subtle prograde travelling pattern identified above 
in He I is tentatively mimicked in C IV, which supports an inter- 
pretation in terms of photospheric velocity fields; though con- 
firmation of this result clearly requires much higher signal-to- 
noise time-series data. 



us I - Jam. 25, 3t)03 C Df - JEOl, 3B, 3003 




Fig. 5. An example of absorption line profile changes in 
Hel i5876and C IV A5801 over ~ 5.3 hours. Individual quo- 
tients with respect the mean are plotted, for the same dynamical 
scale as in Fig. 4. The tick-marks along the right-hand side in- 
dicate the timings of the individual spectra, for both panels. 
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3.3.2. The 2004 time series 

Interestingly, the 0.36-d period derived from the 2002 data 
of a Cam is not present in the observations secured during 
November 2004. The line profile changes are generally weaker 
in 2004 and the Fourier analysis does not reveal any significant 
frequencies in Hel /15876 , Ha or Hel /16678 . However, whilst 
the time-series of 26 and 28 November, 2004 are almost steady, 
a different episodic property is evident during 27 November, 
2004. In contrast to the weak fluctuations within the absorption 
trough seen in 2002 (Figs. 4 and 5), a systematic velocity shift 
occurs in the absorption components of the He I lines during 
this night; this behaviour is shown in Fig. 6. 

The central velocity (and wings) of Hel /16678 shifts by 
up to ~ 30 km s ' during the night, while the motion of 
Hel /I5876has apeak shift of ~ 15 km s ' (but in a coiTespond- 
ing manner). Interestingly, these changes in absorption velocity 
are accompanied by relatively strong changes in the H„ emis- 
sion flux, which are centred at rest velocity and occur within 
Ve sin(/) (see Fig. 6). The H„ flux is a maximum in this region 
when the He I absorption lines have their greatest blue-ward 
motion. 

We stress that the behaviour documented in Fig. 6 is iso- 
lated to November 27, 2004, it was not witnessed for example 
during 2002 or during the remaining nights in 2004. This be- 
haviour (also shown here in Fig. 1) is empirically very differ- 
ent from the broader changes due to stellar wind variability that 
grossly afifect the blue and red wings of Hq. and might in prin- 
ciple be due either to variations in radial velocity of the under- 
lying photospheric profiles or to localised changes in density 
stratification in the innermost part of the wind. 

To test the first possibilit y we emp l oyed t he approx- 
imate method deve loped by iPuls et al.l (119961) (see also 
iMarkova et al l l2004 and applications in Sect. 4) to calculate a 
pair of 'representative' H„ model profiles: One with a normal 
TLUSTY (reff= 30000K, logg= 3.2) photospheric profile and 
the other with the same TLUSTY input but now shifted in ra- 
dial velocity by -30 km s ' (cf. Hel /I6678in Fig. 6). The cor- 
responding model Hq. profiles are shown in Fig. 7, demonstrat- 
ing some correspondence to the observed Hq. changes shown 
in Fig. 6. The effect of shifting the core Hq photospheric ab- 
sorption blue-wards is (naturally) to reveal more of the central 
Hq. emission. 

It seems therefore possible that the photospheric lines of a 
Cam are affected by large-amplitude radial velocity shifts. It 
is unknown whether this characteristic is periodic over a much 
longer time-scale, but the phenomenon was only seen once in 
our 2002 and 2004 data sets. The radial velocity motion is dif- 
ferential, such that the extent of the velocity shifts is differ- 
ent for lines formed at varying depths in the atmosphere. This 
likely reflects a velocity gradient in the atmosphere, though the 
physical origin of this deep-seated disturbance is not known. 

On the other hand, the Hq, variability pattern illustrated in 
the lower right panel of Fig. 1 seems to be quite similar to the 
ones obser ved bv Markova on Dec. 31, 1998 and on Jan 6, 1999 
dMarkova' 2002, see Fig.4 ). In addition, the radial velocity be- 
haviour of the Hel /I66781ine during those two nights in 1998 
and 1999 was similar to that shown in Fig 6, i.e. the increase 



in the Hq emission flux localised within + sin(/) was accom- 
panied by a shift to the blue in the radial velocity of the line. 
The time resolution of the 1998-1999 data set was very low (3 
spectra per night only) and thus prevents a more detailed com- 
parison of the phenomena to the 2004 event. Nevertheless, the 
empirical similarity of the Nov. 27, 2004 variations and those 
observed on Dec. 31, 1989 supports a common physical origin, 
perhaps in terms of large-scal e perturbations in the inner-most 
wind regions, as suggested bv lMarkoval ( l2002h . 

4. Stellar wind characteristics 

It is well known that the stellar winds of OB stars are highly 
variable. Evidence for wind structure is seen in modulated 



UV and optical lines (e.g. 


'Kaperetall Il996t Ide Jong etall 


2001; Prinja etal. 2004; 


Markova et al 




1, excess flux at 



infrar ed and millimetre wavelengths fRunagresj^^lommgJ, 
[1996'), and detection of X-rays (e.g. Chlebowski et al Il989t 
ICassinelli et al, 2001). a Cam is no exception in this respect as 
already demonstrated in Figs. 1 and 2 (see also references in 
Sect. 1). 

The Hel /I58761ine of a Cam offers a weak diagnostic of 
variable wind conditions via fluctuations in the blue-ward ab- 
sorption wing (Fig. 1); unfortunately we only have coiTespond- 
ing Hffdata for the November 2004 run. Excluding the excep- 
tional behaviour seen on 27 November, 2004 (Sect. 3 . 1 ), a max- 
imum increase in Hq equivalent width between -100 to -600 
km s^'of ~ 0.4A is accompanied by a 0.05 A absorption en- 
hancement in He I. We do not find any evidence for periodic or 
modulated stellar wind changes over ~ hourly time-scales. Our 
data set is not suitable for an assessment of temporal behaviour 
on stellar rotational time-scales. 

We focus in this section instead on the large-scale, epoch- 
to-epoch changes evident in Hq, which are, of course, superim- 
posed on the more subtle photospheric and deep atmospheric 
structure described in Sect. 3. Our goal here is to explore ex- 
treme Hq morphologies in a Cam that provide challenges for 
line-synthesis modelling and likely reflect variable density dis- 
tributions in the inner- wind regions. 

The overall HQline profile shape can vary substantially over 
yearly time-scales. A particular example is shown in Fig. 8 
(upper panel) where we compare line profiles from November 
2004 and January 1999. The ~ 30% change in equivalent width 
is mostly around line centre and toward the red-ward emission 
wing. To fit the tw o profil es we employed again the approxi- 
mate method of Pui s et al ](|I996), modified by Markova e t alJ 
(2004), to account for the effects of line-blocking and blanket- 
ing. We find that the observed red-wing and peak are reason- 
ably well matched for fixed terminal velocity = 1550 km s"', 
/3 velocity law index - 1.05 and varying the mass-loss rate (M) 
between 6.5 x lO"*" Moyr^'to 5.1 x lO"''M0yr"'. The model 
with the lowered emission flux cannot however reproduce the 
blue-ward regions of the line profile from 1999, since it under- 
predicts the emission in this region. 

The lower panel in Fig. 8 shows a markedly different tem- 
poral characteristic in the Hq profile of a Cam. These two spec- 
tra were taken 4 days apart in January 1999, and they exhibit 
broadly similar red-ward emission and peak fluxes. However, 




Fig. 6. Upper panels: Montage of quotient spectra with respect to the mean, showing variabiHty in Hq., Hel /I6678and 
Hel /15876 during 27 November 2004. Lower panels: Representative examples of the original (not residual) line profiles. Note the 
radial velocity shifts in the He lines and the corresponding change in H„peak emission flux. (For Hel /15876 the Na I D interstellar 
lines (Sect. 2) provide an excellent wavelength fiducial.) 



in this case the blue-ward regions of the profile are highly vari- 
able between ~ -50 to -450 km s"' (i.e. well beyond the do- 
main of He II /16560.2). Whilst models with M differing be- 
tween -6.2 X lO-*" Moyr-i and 5.9 x lO-*" Moyr-' (v^ = 1550, 
j6 = 1 .05) can broadly match the red-wing, they once again can- 
not reproduce the dififering morphologies of the blue-ward re- 
gions. 

We note therefore that though the variable red-ward 
H„ wings allow some constraint on M and /3, the blue-ward 
fluctuations pose a more serious discrepancy for spherically 
symmetric stellar wind models with smooth density stratifica- 
tion. 



i) systematic changes in the photospheric lines, with some ev- 
idence that the behaviour may be linked to surface velocity 
fields due to non-radial pulsations, 

ii) short time-scale (~ hourly) perturbations that differentially 
shift the radial velocity of He I lines, and potentially give 
rise to localised changes in Haemission, 

iii) a variable stellar wind, which may be parameterised typi- 
cally in terms of a ~ 5% change in mass-loss rate over sev- 
eral days, but very occasionally this difference can rise to 
~ 30%. The red-ward emission changes are accompanied 
by more enigmatic blue-ward changes that may result from 
absorption effects linked to localised density structures in 
the outflow. 



5. Discussion 

We have demonstrated in this study that the primarily, or par- 
tially, wind-formed spectral lines of the late O supergiant a 
Cam can exhibit variability that arises in different physical 
regimes of the star. This 'contamination' of wind signatures 
by variable non-wind components is directly relevant to line- 
synthesis modelling, particularly in the context of revising 
mass-loss rates due to the presence of small- and large-scale 
structure in the outflows. In particular, we have found that the 
temporal behaviour of a Cam's optical spectrum is charac- 
terised by at least three different variability patters: 



Hot massive stars are expected to be subject to photospheric 
instabilities and/or pulsations. During the last decade the search 
for non-radial pulsations in OB stars has been largely moti- 
vated by the need to test the hypothesis of the "photospheric 
connection" which aims at linking the cyclic variability ob- 
served in the winds of these stars to processes in the stellar 
photosphere. The task however turned out to be observation- 
ally very demanding and as a result only a limited number of 
OB stars ar e currently flagged as confirmed or suspected pul- 
sators re.g. iFullerton. Gies fc RoltonL jl 99(4 iReid & Howartll 
fl99aide Jong et^EoOltlKauferltalEoO^ 

The results derived in this study raise the possibility that the 
supergiant a Cam might also be a non-radial pulsator. Radial 
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Fig. 7. 'Representative' model Ho-profiles, demonstrating the 
effect of shifting the TLUSTY photospheric models by -30 km 
s"' (dashed line; see Sect. 3.3). 
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Fig. 8. H„model fits (solid and dashed lines) to extreme 
changes observed in 2004 and 1999 (filled and open circles). 

velocity variations in optical absorption Unes of a Cam have 
been rep orted previ ouslv bv IZeinalov & MusaevI il986 ) and 
iMusaev & Snezhkol( ll988l) who argued that these variations are 
likely periodic with a possible cycle length of 1.35, 1.54 or 
2.93 days (with no preference to any of them). However, our 
analysis of the 2002 and 2004 time-series did not support any 
of these periods but instead revealed a new one of 0.36 days, 
which may be interpreted in terms of prograde low-order non- 
radial pulsational behaviour. The ~ 8-9 hour modulation may 



persist over two months in 2002, albeit with different phase 
shifts. This signal is not evident in data collected more than 
2 years later (November 2004); it is also not seen in the stel- 
lar wind components of the line profiles. The interpretation of 
the 0.36-day variation in terms of photospheric velocity fields 
is also supported by weak fluctuations in the shallow photo- 
spheric lines of C IV AA 5801, 5812. 

Apart from the 0.36-day periodic variation detected be- 
tween Jan- March 2002, our observations also documented oc- 
casional large-amplitude (of about 15 to 30 km s"') differential 
radial velocity shifts in Hel /15876 and Hel /16678 absorption 
lines, accompanied by an increase in the H^emission flux cen- 
tred at Vsys and localised within + v^. sin(/). 

Large-a mplitude differentia l velocity shifts have been dis- 
covered bv lPrinia et all ( l2004 in the optical absorption lines 
of e Ori (BO la) and they are also present in the photospheric 
lines of /3 Cepheid stars. However in both these cases the ve- 
locity shifts are regular and periodic while those seen in a Cam 
occurred over at least several hours, during a specific night in 
2004, with no recurrence during the other nights of the run. 
This phenomenon can be interpreted either in terms of large- 
amplitude changes in photospheric velocity fields or as the re- 
sult of deep-seated, large-scale wind structure that occasionally 
modulates the density distribution in the innermost part of the 
wind. 

Discrete Absorption Components moving from red to blue 
across the unsaturated UV resonance lines are the characteristic 
signature of large scale structure in O star winds. Unfortunately 
the UV resonance lines of a Cam are saturated and do not 
provide any detai led information on variable wind features 
(lKaperetaill997l) . We have re-examined the lUE time-series 
of a Cam from February 1991 and December 1994. Whilst 
there is subtle evidence for low velocity fluctuations in Si IV 
AA1400 (and Si 111 /11207), the clearest variability is in the soft 
blue wings of the saturated UV lines. For example, during the 
1991 lUE observing run of ~ 5 days, additional absorption 
is evident in the blue wing of Si IV beyond Vm (- 1550 km 
s"') out to ~ -1700 km s"'. This enhancement persists for ~ 
24 hours and does not recur during the following 4 days. This 
finding is at least consistent with the notion based earlier on op- 
tical observations, that large scale structure may develop in the 
wind of a Cam; the present study suggests however that these 
structures are somewhat unstable and thus do not give rise to 
longer-term periodic phenomena in the wind. 

The possibility that the wind of a Cam is not smooth 
but structur ed has got additional support by t he rec ent results 
derived by iFullerto n. Massa & Prinia (2006) and Puis e t all 
(|2006) . In particular, Fullerton. Massa & Prinia C2006) have 
modelled the P V lines of a Cam using the SEI method to derive 
products of mass-loss rate and ionization fraction. Adopting 
mass-loss rates derived from radio observations, they find re- 
sulting ionization fractions, q(P'*^), = 7x10"^ for a Cam. In the 
case of a Cam the low ion fraction is highly discrepant with re- 
sults from NLTE stellar atmosphere models which predict that 
the P V ion fraction should be dominant and thus close to unity 
at these effective temperatures. The discrepancy is thought to 
indicate that the wind is highly clumped and has very small 
volume filling factors. 
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Further evidence in support of the clumped structure of 
a Cam wind has been derived bv lPuls et alJ (l2006l) . By means 
of a simuhaneous model analysis of H„, IR and radio data these 
authors showed that the wind of a Cam likely has a clumping 
factor in the innermost region ( between 1 .05 ^ and 27?»), being 
a factor of 2.6 larger than in the outermost wind (beyond 50/?,). 
Using the H|,emission wings as a tracer of wind clumping in 
the intermediate wind (between 2 and 15/?,), Puis et al. further- 
more found that the clumping factor in this region is similar or 
fractionally larger than in the inner wind. 

Interestingly, the prop erties of t h e TVS of HQ,(derived in 
the present study and by 'Markova' ("2002")) also seem to be 
consistent with results obtained vJji_line-profile simulations in- 
cluding wind clumps jMarkova et ali 1200^ . In particular, the 
velocity limits as well as the established asymmetry of the 
TVS are both consistent with a model including 'broken shells" 
(i.e. clumps), uniformly distributed over the wind volume, 
with density contrast 6p/po of +0.35 (model series BS21 in 
(Markova et al, 2005)) where po is the stationary density. Thus, 
it may be that Hu-night-to-night variability observed in a Cam is 
at least partially connected with presence of small-scale struc- 
ture (clumps, blobs) in the innermost part of the wind. 

Overall, there is now substantial observational evidence 
which suggests that the wind of a Cam is not smooth but struc- 
tured over, both, large and small spatial scales. In view of this 
possibility it is therefore not surprising that we failed to repro- 
duce the observed set of H„profiles (Sec. 4) by means of model 
calculations. Indeed, the approximate method of Puis et al. 
ill 9961) employed here, as well as most of the state-of-art model 
atmosphere codes, rely on the assumption of a globally station- 
ary wind with a smooth density /velocity stratification, coupled 
to a 'stable' photospheric line contribution, and are therefore 
incapable in principle of describing the observed blue-ward 
(absorptive) behaviour of H^in a Cam and of any other star 
with a substantially structured wind. 

Acknowledgements. This study was in part supported by funds 
from the National Scientific Foundation to the Bulgarian Ministry 
of Education and Science (F- 1407/2004). We are also grateful to 
Bulgarian Academy of Sciences and the Royal Society (UK) for a 
collaborative research grant. NM acknowledges the financial support 
of the Bulgarian Academy of Sciences and PPARC. NM is also grate- 
ful to her colleague Sabotinov (BAS) for his understanding and moral 
support during this study. 



Fullerton, A. W., Gies, D.R & Bolton, C.T. 1996, ApJS 103, 475 
Fullerton, A.W., Massa, D., & Prinja, R.K. 2006, ApJ, 637, 1025 
Hamann, W.-R. 1980, A&A 84, 342 
Haser, S.M. 1995, PhD thesis, Univ Munich 
Heger A., et al. 2003, ApJ, 591, 288 
HiUier, D.J., Miller, D.L., 1998, ApJ, 496 

Hoffleit, D. & Jaschek, C. 1982, The Bright Star Catalogue, New 

Haven: Yale University Observatory (4th edition) 
Kaper, L., Henrichs, H.F., Nichols, J.S., Snoek, L.C., Volten, H., & 

Zwarthoed, G.A.A. 1996, A&AS, 116, 257 
Kaper, L., et al. 1997, A&A, 327, 281 

Kaufer, A., Stahl, O., Prinja, R.K., & Witherick, D., 2006, Ac&A, 447, 
325 

Lamers, H.J.G.L.M., Snow, Th., de Jager, C, & Langerwere, A. 1988, 

ApJ, 325, 342 
Markova, N., & Valchev, T. 2000, A&A, 363, 995 
Markova, N. 2002, A&A, 385, 479 

Markova, N., Puis, J., Repolust, T, & Markov, H. 2004, A&A, 413, 
693 

Markova, N., Puis, J., Scuderi, S., & Markov, H. 2005, A&A, 440, 
1133 

Morel, T, Marchenko, S.V., Pati, A.K., Kuppuswamy, K., Carini, 
M.T, Wood, E., Zimmerman, R. 2004, A&A, 215, 45 

Musaev, EA. & Snezhko, L.I. 1988 SvAL, 14, 68 

Pauldrach, A.W.A., Hoffmann, T.L., & Lennon, M. 2001, A&A, 375, 
161 

Penny, L.R. 1996, ApJ, 463, 737 

Prinja, R.K., Massa, D., & Fullerton, A.W. 2002, A&A, 388, 587 
Prinja, R.K., Rivinius, Th., Stahl, O., Kaufer, A., Foing, B. H., Cami, 

J., & Orlando, S. 2004, A&A, 418, 727 
Puis, J. et al. 1996, A&A, 305, 171 

Puis, J., Urbaneja, M. A., Venero, R., Repolust, T, Springmann, U., 

Jokuthy, A., & Mokiem, M. R. 2005, A&A, 435, 669 
Puis, J., Markova, N., Scuderi, S., et al., 2006, A&A, submitted 
Reid, A.H.N., & Howarth, I.D. 1996, A&A, 311, 616 
Repolust, T, Puis, J., & Herrero, A., 2004, A&A, 415, 349 
Roberts, D.H., Lehar, J., & Dreher, J.W. 1987, AJ, 93, 968 
Runacres, M.C., & Blomme, R. 1996, A&A, 309, 544 
Schaller, G., Schaerer, D., Meynet, G., & Maeder, A. 1992, A&AS, 
96, 269 

Zeinalov,S.K & Musaev, EA. 1986 SvAl 12, 125 
Walborn, N.R. 1973, AJ, 78, 1067 



References 

Bouret, J.C., Lanz, T, & HiUier, D.J. 2005, A&A, 438, 301 
Cassinelli, J.P, Miller, N.A., Waldron, W.L., MacFarlane, J.J., & 

Cohen, D.H. 2001, ApJ, 554, L55 
Chlebowski, T, Harnden, ER.J., & Sciortino, S. 1989, ApJ, 341, 427 
de Jong, J.A., et al. 2001, A&A, 368, 601 
Ebbets, D. 1982, ApJS 48, 399 
Fullerton, A.W. 1990, Ph.D. thesis, Univ. Toronto 

^ The wind below 1.05/?, was assumed to be unclumped in agree- 
ment with results from hydrodynamical simulations which showed 
that the line-driven instability needs some time to grow before sig- 
nificant structure is formed. 



